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ABSTRACT 



This paper extends an earlier high-temperature lattice calculation of the renormalized 



^ ■ Green's functions of a D-dimensional Euclidean scalar quantum field theory in the Ising 
limit. The previous calculation included all graphs through sixth order. Here, we present 
the results of an eleventh-order calculation. The extrapolation to the continuum limit in 
the previous calculation was rather clumsy and did not appear to converge when D > 2. 
Here, we present an improved extrapolation which gives uniformly good results for all real 
values of the dimension between D = and D = 4. We find that the four-point Green's 
function has the value 0.62 ± 0.007 when D = 2 and 0.98 ± 0.01 when D = 3 and that 
the six-point Green's function has the value 0.96 ± 0.03 when D = 2 and 1.2 ± 0.2 when 
D = 3. 



PACS numbers: ll.10.Ef, ll.10.Jj, ll.10.Np 
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There have been many attempts to calculate the coefficients in the effective potential 
of a Euclidean scalar quantum field theory in the Ising limit. These coefficients are just the 
dimensionless renormalized 2n-point Green's functions evaluated at zero momentum. Tech- 
niques that have been used to determine these Green's functions include high-temperature 
lattice expansions, Monte Carlo methods, and epsilon expansions. For the case D = 3 a 
complete list of references is given in a recent paper by Tsypin. 1 

In a series of papers 2 ' 3 high-temperature lattice techniques were used to obtain the 
dependence of the Green's functions upon the Euclidean space-time dimension D for D 
ranging continuously between and 4. In this work a strong-coupling calculation to sixth 
order was performed analytically on a hypercubic lattice in D dimensions and Pade extrap- 
olation techniques were invented to obtain the continuum limit. 4 The current paper is an 
addendum to the work in Refs. 2 and 3. Here, we extend the strong-coupling calculation to 
eleventh order. Furthermore, we use an improved Pade extrapolation method that relies 
on information taken from the results of a large- N calculation and our recent studies of 
dimensional expansions for quantum field theory. 5 ' 6 ' 7 

Our strong-coupling lattice calculations are identical to those described in Ref. 2 ex- 
cept that the graphs were generated using a FORTRAN program and evaluated analytically 
using MACSYMA. The eleventh-order calculation involves several hundred times as many 
graphs as the sixth-order calculation. We have verified the accuracy of our expansions 
for the specific cases of D = 2, D = 3, and D = 4 dimensions by comparing them with 
previous calculations 8 . 

The lattice series for the renormalized four-point and six-point Green's functions are 
as follows: 
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y-D/2 

74 = — — [ 1 + 4 D y + (4 D 2 - 10 D) y 2 + 16 D y 3 + (30 D - 80 L> 2 ) y 4 

+ (256 L> 3 + 104 D 2 - 192 D) y 5 

+ (-704 D 4 - 1736 D 3 + 2508 D 2 - 656 D) y 6 

+ (1792 D 5 + 10432 D 4 - 11232 D 3 - 3872 D 2 + 4992 D) y 7 

+ (-4352 D 6 - 45600 D 5 + 11456 L> 4 + 123672 D 3 

- 128440 D 2 + 35542 D) y 8 

' (1) 

+ (10240 D 7 + 168320 D 6 + 181248 D 5 - 1052576 D A 
+ 2615584/3 D 3 + 76664 D 2 - 681472/3 D) y 9 

+ (-23552 D 8 - 558208 D 7 - 1630272 D 6 + 5391904 D 5 

- 1011536/3 D 4 - 10102936 D 3 + 29622092/3 D 2 - 2720752 D) y 10 

+ (53248 D 9 + 1718272 D 8 + 9081856 D 7 - 18274816 D 6 

- 113682176/3 D 5 + 367432576/3 D 4 - 292976128/3 D 3 
+ 18425408/3 D 2 + 14757984 D) y 11 + . . . ] 

and 

y-D 

7 6 = [ 1 + 6Dy + (12 D 2 - 6D)y 2 + (8 D 3 - 12 D 2 - 20 D) y 3 

30 

+ (48 D 2 + 48 D) y A + (-96 D 3 - 816 D 2 + 528 D) y 5 

+ (192 D 4 + 4640 D 3 - 2736 D 2 - 560 D) y 6 

+ (-384 D 5 - 18432 D A - 10800 D 3 + 46512 D 2 - 23040 D) y 7 

+ (768 L> 6 + 61440 D 5 + 188352 D 4 - 510816 L> 3 + 357324 D 2 - 72492 D) y 8 

+ (—1536 D 7 - 184576 D 6 - 1274880 D 5 + 2653440 D 4 

- 77496 D 3 - 2911496 D 2 + 1698240 D) y 9 

+ (3072 D 8 + 517632 D 7 + 6280704 D 6 - 6584832 D 5 - 27745840 D 4 
+ 65401176 D 3 - 49332608 D 2 + 11853912 D) y 10 

+ (-6144 D 9 - 1382400 D 8 - 25928448 D 7 - 13343232 D 6 + 286690784 D 5 

- 516057392 D 4 + 211594432 D 3 + 210150872 D 2 - 153291336 D)y 11 + ...], 

(2) 
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where y = (Ma) -2 , a is the lattice spacing, and M is the renormalized mass, which is 
obtained from the two-point function as explained in Ref. 2. 

The quantity y/y is the dimensionless correlation length. The continuum limit a — > 
corresponds to infinite correlation length. To obtain the continuum Green's functions it 
is necessary to extrapolate the formulas in (1) and (2) to their values at y = oo. Direct 
extrapolation to the continuum limit of either series in (1) or (2) leads to a sequence 
of extrapolants that becomes badly behaved when D increases beyond 2; we find that 
extrapolations as functions of D do not converge to a limiting curve (see Figs. 1 and 2). 
However, for D near these extrapolants are well behaved and converge rapidly to the 
known exact values 2 74 = 1/4 and 76 = 1/4 at D = 1. 

To improve our extrapolation we make the following observation. We consider a scalar 
field theory having an O(N) symmetry. The model we have studied above corresponds to 
the case N = 1. In the limit iV — > 00 one can solve for the Green's functions exactly. We 
obtain the following lattice results: 



N (n=oo) = y D/2 rri 

74 4/ oo dtte-*[e- 2 ^/ (2t?/)P 1 ' 

and 

^2 (;v=oo) = y~ D ^dtt\-t[e-^I G {2ty)] D 

76 l2 (/ 00 rf^^e- ^ [e- 2t no(2^ 2/ )]^) 3, 

where we have summed over the external indices. In the continuum limit y — > 00 we have 

74 4 r(2-D/2) {) 



and 



N i <*=«,) _ (4 7 r)^r(3 - g/2) 

76 " 12 [r(2-D/2)] 3 ' 1 J 



where D lies in the range < D < 4. Each of these functions rises from its value at D = 0, 
attains a maximum, and falls to at D = 4. 

Under the assumption that the Green's functions for N = 1 vanish at D = 4 like those 
in (5) and (6) we can extract such a behavior from the series (1) and (2) by performing a 
Borel summation as follows. Consider the lattice series in (1) for 74. This series has the 
general form 

— D /2 00 

74 = (7) 

where Pk(D) are polynomials of maximum degree k. One can read off the first eleven 
polynomials Pk(D) from (1). We can rewrite (7) as 

where Qk(D) is another polynomial in D. Next, we insert the identity 

1 r°° 

1 = (F+T)! / * tt+lf "' < 9 > 

for each term in the sum. This converts (8) to the form 

1 r°° 

— = dtf-We-tfiyt), (10) 
1^74 Jo 



where we define 



m=x SifcW 



fc=0 



(11) 



fc=0 

where again Rk(D) are polynomials in D. 

We now take the continuum limit of the expression (10). Assuming that /(oo) exists 
in the limit y — > 00 so that (10) separates into a product of two terms, we have 



Thus, we have forced the continuum limit of the four-point function to take the form of 
the large- N result in (5) apart from /(oo), which is a smoothly varying function of D. 

There is an immediate indication that the Borel summation leading to (12) has a 
significant impact. We find that the polynomials Rk(D) in (11) are significantly simpler 
than the original polynomials Pk(D) in (7); the polynomials Rk(D) have maximum degree 
[k/2] — 1, about half the degree of the polynomials Pk(D). The polynomials Rk(D), 
k = 0, 1, 2, . . . , 11, are as follows: 

Ro(D) = 1, 
Ri(D) = —2, 
R 2 (D) = 11/3, 
R 3 (D) = -16/3, 

„ 73 „i 1346 „ 76691 

Be(0) = 2838 B " W°+ 11340' < 13 > 

7V ' 1134 1260 2835 

Rs(D) = _ D 3 + JW^L _ ?837 2627137 
8V y 8100 170100 3240 544320 

d /rj\ _ 163 n3 14699 2 17111 235681 
^ 9(jD) -3150 D "14175^ + W D ""64800"' 

„ 457 . 1338257 o 2015579 ~ 9257497^ 230357209 

R W (D) = D D 3 H £) D H , 

1UV ; 267300 5613300 1403325 3207600 89812800 ' 

„ 611 . 6066953 ~ 3923681 ~ 110092603^ 75132389 

Rn(D) = D 4 H D 6 D H £) . 

V ; 69300 11226600 2245320 44906400 44906400 



The resummation of the lattice series as performed above reduces the problem of 
extracting the continuum limit to finding the value of /(oo). This is done using the same 

6 



Pade techniques as were used in Ref. 2. If we perform this numerical calculation we obtain 
a sequence of approximants, one for each new order in perturbation theory. The first eleven 
such approximants for 74 in (12) are plotted in Fig. 3. Each approximant is a continuous 
function of D for < D < 4. Note that the approximants are smooth and well behaved; 
the sequence is monotone increasing and appears to converge uniformly to a limiting curve. 
The dotted line on Fig. 3 is an extrapolation of these approximants to this limiting curve 
obtained using Richardson extrapolation. 9 

There are a number of ways to assess the accuracy of the limiting curve. First, one 
can Taylor expand this limiting curve about D = as a series in powers of D. This Taylor 
series has the form 

^imiting curve ^ = jL ( 1 + 1 £ + Q.64 I) 2 + 0. 19 L> 3 + 0.03 L> 4 + . . .) . (14) 

We may then compare this Taylor series with that recently obtained 7 using dimensional 
expansion methods: 

^imensional expansion^ = _L(1 + \.\%D + 0.62 D 2 + 0.18 D 3 + 0.03 D 4 + ...). (15) 

J. — ■ 

Note that the coefficients of these two series are almost identical. Second, we can examine 
the limiting curve at D = 1, for which the exact value 74 = 1/4 is known. At this value of 
D the limiting curve has the value 0.2526 so it is slightly high by about 1%. 

In Figs. 4 and 5 we demonstrate how we obtain the limiting curve for the cases D = 2 
and D = 3. We have plotted the nth-order Richardson extrapolants for the approximants in 
Fig. 3 for n = 1, 2, . . . , 5 versus the inverse order of the approximants. We then determine 
where each of these extrapolants crosses the vertical axis (each intersection is indicated by 
a horizontal bar). Finally, we extrapolate to the limiting value of these intersection points. 
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This procedure gives the value 74 = 0.620 ± 0.007 at D = 2, indicated in Fig. 4 by a fancy 
square. This result is to be compared with 74 = 0.6108 ± 0.0025 obtained by Baker and 
Kincaid 10 . Similarly, in Fig. 5 we find that the limiting curve gives 74 = 0.986 ± 0.010 
at D = 3. This value compares reasonably well with previous results, as tabulated in 
Ref. 1. For example, Baker and Kincaid 11 obtain 0.98, Monte Carlo studies give results 
between 0.9 and 1, and renormalization group studies give results around 0.98. Note that 
the limiting curve in Fig. 3 has a maximum extremely close to D = 3; numerically, the 
maximum occurs at D = 3.03. 

The same procedure that was used to extrapolate (1) to the continuum and thereby to 
obtain a plot of 74 as a function of D can be applied to (2). We perform a Borel summation 
of the series in (2) as follows. The lattice series in (2) for 7 6 has the general form 

-d 00 

7 6 = ^E^(^' ( 16 ) 

fc=0 

where Sk(D) are polynomials of maximum degree k. One can read off the first eleven 
polynomials Sk(D) from (2). From the structure of 7g Ar_0 °' > in (4) we are motivated to 
rewrite (16) in the form 

y D / 2 j oo dtt 2 e- t [e- 2t yi (2ty)] 
30 76 

where Tk(D) is another polynomial in D. Again, we insert the identity (9) for each term 
in the sum in (17). This converts (17) to the form 



k=0 



(17) 



y D li^dtee-\e-^h{2ty)\ D 
30 76 



" />oo "1 3 

Jo 



(18) 



fc=0 



where we define 

r. k 

^ (k + i)T 

,D/2 (19) 

xJ2Uk(D)x k ) , 

A:=0 / 

where Uk(D) are polynomials in D of degree [k/2] — 1 similar in structure to those in (13). 

Next, we take the continuum limit of the expression (18). Assuming that g(oo) exists 
in the limit y — > oo we find that (18) separates into a product of several terms and we have 

= (4^-^(3 -D/2) 
16 30[r(2-D/2)(7(oo)] 3 ' 1 ' 

Thus, we have forced the continuum limit of the six-point function to take the form of the 
large- N result given in (6) apart from g(oo), which is a function of D. 

Again, the resummation of the lattice series reduces the problem of extracting the 
continuum limit to finding the value of g(oo). This is done using the same Pade techniques 
as were used in Ref. 2. We perform this numerical calculation and obtain a sequence 
of approximants, one for each new order in perturbation theory. The first eleven such 
approximants for 7 6 in (20) are plotted in Fig. 6. Each approximant is a continuous 
function of D for < D < 4. As in Fig. 3 the approximants are smooth and well behaved; 
the sequence is monotone increasing and appears to converge uniformly to a limiting curve 
indicated in Fig. 6 by a dotted line. This limiting curve is again obtained using fifth- 
order Richardson extrapolation. The limiting curve at D = 1 passes through the value 
76 = 0.240 which differs from the exact value 76 = 1/4 by about 4%. 

The limiting curve predicts that 7 6 = 0.96 ± 0.04 at D = 2 and 7 6 = 1.2 ± 0.1 at 
D = 3. This value is lower than most previous results, as tabulated in Ref. 1, but it is 
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1,12 



certainly larger than zero. By comparison, an epsilon expansion around D = 4 gives 

Tfi 20 o o 

2e - — e 2 + 1.2759 e 3 + . . . . (21) 



( 74 ) 2 27 

This series appears to be divergent but a direct optimal truncation of the series after one 
term with e = 1 gives the value 76 = 1.9 ± 0.7. (Here, we have substituted the value 
74 = 0.986 given above.) However, if we perform a (1, 1)-Pade summation of this series, 
which seems justified because of the alternating sign pattern, we obtain the smaller value 
76 = 1.66 ± 0.28, in better agreement with our predicted value. 

Finally, we observe that the maxima of 72^ as a function of D appear to follow a 
pattern. We observed already that 74 has a maximum that is close to D = 3. Here, 
we find that the limiting curve for 75 has a maximum at D = 2.66 which is very close 
to 8/3. An interesting conjecture is that in general the maximum might be located at 
Dmax = 2 j ^ ) the value of D for which a <^> 2l+2 theory becomes free. 

We thank the U. S. Department of Energy for funding this research. 
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FIGURE CAPTIONS 

Figure 1. The first k Pade approximants, k = 1, 2, 11, to the continuum limit of the 
four-point Green's function. These approximants are constructed from the lattice series in 
(1) using the Pade procedure explained in Ref. 2. Observe that the approximants are well 
behaved when D is near 0; in particular, they converge nicely to the exact value 1/4 at 
D = 1 (indicated by a plus sign). However, when D increases beyond 2 the approximants 
become irregular and do not seem to converge to a limiting function of D. (Some of the 
approximants reach zero and terminate as D increases because they become complex.) 

Figure 2. Same as in Fig. 1 except that we have plotted the first k approximants to the 
six-point Green's function for k = 1, 2, 11. Again, we find that the approximants 
converge to the exact value 1/4 at D = 1 but that they are irregular for D > 2. 

Figure 3. First eleven approximants to 74 plotted as functions of D for < D < 4 [see (12)]. 
The approximants form a monotone increasing sequence of curves; the labeling indicates 
the order. Note that the approximants are smooth curves that seem to be tending uniformly 
to a limiting curve. We have obtained this limiting curve (dotted curve) by means of fifth- 
order Richardson extrapolation 9 . The exact result 74 = 1/4 at D = 1 is indicated by a 
plus sign; the limiting curve passes within 1% of this point. 

Figure 4. Plot of the nth-order Richardson extrapolants for the approximants in Fig. 3 
for n = 1, 2, . . . , 5 versus the inverse order of the approximants for the case D = 2. 
Each Richardson extrapolant is linearly extended (dash-dotted line) until it intersects the 
vertical axis. Each intersection is indicated by a horizontal bar. We then extrapolate to 
the limiting value of these intersection points, indicated by a fancy square. This procedure 
predicts that 74 = 0.620 at D = 2. 
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Figure 5. Same as in Fig. 4 except that D = 3. This extrapolation procedure predicts that 
74 = 0.986 at D = 3. 

Figure 6. First eleven approximants to 76 plotted as functions of D for < D < 4 [see 
(20)]. The approximants form a monotone increasing sequence of curve as indicated by 
the labeling. As in Fig. 3 for the four-point function, the approximants are smooth curves 
that seem to be tending uniformly to a limiting curve. This limiting curve (dotted curve) 
is a fifth-order Richardson extrapolation. The exact result 76 = 1/4 at D = 1 is indicated 
by a plus sign; the limiting curve passes within 4% of this point. 
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